Electronic structure and orbital polarization of LaNiOs with a reduced coordination 

and under strain: first-principles study 
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First-principles density functional theory calculations have been performed to understand the 
electronic structure and orbital polarization of LaNiOa with a reduced coordination and under strain. 
From the slab calculation to simulate [001] surface, it is found that d^,,2_r2 orbital occupation is 
significantly enhanced relative to d^2_y2 occupation owing to the reduced coordination along the 
perpendicular direction to the sample plane. Furthermore, the sign of the orbital polarization 
does not change under external strain. The results are discussed in comparison to the bulk and 
heterostructure cases, which sheds new light on the understanding of the available experimental 
data. 

PACS numbers: 73.20.-r, 75.70.-i, 71.15.Mb 
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I. INTRODUCTION 

Understanding transition metal oxides is of perpet- 
ual interest and importance in condensed matter physics 
and material science due to their great scientific and 
technological potential [ij . Recent advances in layer-by- 
layer growth techniques of these compounds in the form 
of thin film and heterostructures have created consid- 
erable interest and enabled a further understanding of 
these systems 0]. Exotic material phenomena that are 
clearly distinctive from the 'normal' phases include in- 
terface superconductivity , charge [j, Q and orbital re- 
construction 0, and room-temperature ferromagnetism 
■ Strain-mediated [8l| and electric field control y\ of the 
metal-insulator phase transition are the other examples 
of thin-film phenomena. 

One of the most interesting classes of transition-metal 
oxide systems is thin film LaNiOs (LNO) and het- 
erostructures of LNO and wide-gap insulators such as 
LaAlOs (LAO). Bulk LNO is the only member in the 
i^eNiOa series {Re: La, Pr, Nd,...ete.) that remains 
metallic at all temperatures [l|, As a result, thin 

film LNO and LNO-based heterostructures have been 
studied under various strained conditions to control the 
correlation strength and the metal-insulator phase tran- 
sition [TTl - [l3 | in an effort to find possible high temper- 
ature superconductivity .15.] and conductivity enhance- 
ment [16|. Strain affects most strongly the orbital char- 
acter of the nickelates Recently, several pa- 
pers have focused on how changes in orbital polarization 
in heterostructures and thin films affect the electronic 
structure and thereby the macroscopic material proper- 
ties SfBil. 
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An important issue is the relationship between strain 
and the orbital polarization in various structural classes 
of nickelates. For low-spin trivalent nickelates, the orbital 
polarization P can be defined simply as the difference of 
the two Cg orbital occupations (without renormalization 



P 



(1) 



By definition, a positive orbital polarization indicates a 
higher electron occupation in the dx^-y2 orbital com- 
pared to the d^z2_j.2 state. Recent x-ray linear dichro- 
ism (XLD) measurements on nickelate superlattices and 
thin films have reported intriguing behavior of the or- 
bital polarization P under strain. Beside the transfor- 
mation into a different structural phase [TtI], the polar- 
ization change under compressive and tensile strain has 
mainly been analyzed using Madelung energies. How- 
ever, the detailed changes in electronic structure cannot 
be understood from such a simple electrostatic approach 
fvA [isj . For example, both in film and in heterostruc- 
tures of LNO, compressive strain make the (i3^2_r2 or- 
bital more occupied, as one would expect. However, ten- 
sile strain does not make the polarization go the other 
way which is in a clear contradiction to Madelung pic- 
ture. 

In this paper, we investigate the electronic structure 
and orbital occupation by using first-principles density 
functional theory (DFT) calculations. Special attention 
has been paid to the surface Ni state as its contribu- 
tion might be significant in thin film and in surface- 
sensitive measurements. The electronic structure and 
orbital property are compared to the results of het- 
erostructure and bulk counterpart. Translational sym- 
metry breaking along the direction perpendicular to the 
sample and the reduced coordination at the surface are 
found to strongly change the orbital polarization with re- 
spect to the bulk. The effects of strain are also studied. 
Our results shed new light on understanding the available 
experimental data. 
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FIG. 1: (Color online) Unitcell structure of the double layer 
slab calculation. Large (red), middle (gray), and small (blue) 
size spheres represent La, Ni, and O atoms, respectively. Ni(l) 
and Ni(2) corresponds to the surface-like and bulk-like Ni sites 
respectively. 




o 

Q 



1.2 


.(<=) 






0.8 










0.4 




. 1,. 


r 




n 


A' , 







-6 -4 -2 2 4 

Energy (eV) 



S 1 -4 
in 

1.0 
0.6 



0.6 







1 -0.5 


0.5 1 











-1 -0.5 0.5 1 
Energy (eV) 



FIG. 2: (Color online) (a) Projected DOS and (b) integrated 
DOS of Ni Cg states in (LNO)i/(LAO)i heterostructure. (c) 
Projected DOS and (d) integrated DOS of Ni eg states in 
single layer LNO slab. Red (solid) and blue (dashed) lines 
correspond to d^^i_^i and d^-2_yi states, respectively. The 
Fermi level is set to be zero. 



II. COMPUTATIONAL DETAILS 

For the band-structure calculations, we employed 
TrouUier-Martins type norm-conserving pseudopo- 
tential [23| with a partial core correction and 
linear-combination-of-localized-pseudo-atomic orbitals 
(LCPAO) [l^] as a basis set. In this pseudo-potential 
generation, the semi-core 3p electrons for transition 
metal atoms were included as valence electrons in order 
to take into account the contribution of the semi-core 
states to the electronic structures. Basis orbitals are 
generated by a confinement potential scheme [24] with 
the cutoff radius 7.0 a.u., 5.5 a.u., and 5.0 a.u. for La, 
Ni, and O, respectively. We adopted the local density 
approximation (LDA) for exchange-correlation energy 
functional as parametrized by Perdew and Zunger 
(25j . and used energy cutoff of 400 Ry and k-grid of 
10 X 10 X 10. All the DFT calculations were performed 
using the DFT code OpenMX [H. 

To simulate the [001] surface, the single and double lay- 
ered slab geometries have been investigated. The unit- 
cell structure for the double layer slab is presented in 
Fig. [TJ This is thinner than what is commonly used 
experimentally, but it allows us to obtain a better un- 
derstanding of the surface effects. Vacuum layer thick- 
ness is of 7.0 A which is large enough to make the basis 
function overlap negligible. While the polar surface may 
not be well stabilized in general, these instability can 
be compensated by the large covalency in the nickelates 
p,l, 12, 17, 18]. The geometry relaxation has been per- 
formed with the force criterion of lO^'^ Hartree/Bohr. 
During the relaxation process, the in-plane lattice con- 
stant is fixed considering the substrate effect in the ex- 
perimental situation. The presence of vacuum layer natu- 
rally allows the atomic movement along the out-of-plane 



axis and the volume change. We considered the para- 
magnetic phase because the bulk LNO is paramagnetic 
metal and there is no report on the magnetic ordering 
in the thin film LNO. The orbital polarization P can be 
calculated by integrating the projected density-of-states 
(DOS) up to Fermi level. For the atomic charge den- 
sity we used MuUiken analysis based on the generated 
pseodo-atomic basis orbitals [l^, • 



III. RESULT AND DISCUSSION 

In an ideal cubic structure of LNO, Ni'^^ ion has (for- 
mally) singly occupied degenerate Cg orbitals around the 
Fermi level. This degeneracy can be lifted by, for ex- 
ample, making a heterostructure with a LAO interlayer 
[l5|. It has also been noted that the Cg orbitals in 
LNO/LXO (X=:B, Al, Ga, In) superlattice are positively 
polarized without strain [19|. The positive polarization 
is attributed to the reduced hybridization along the z- 
direction (perpendicular to the sample plane) due to the 
presence of the LAO layers. As a result, the width of the 
bands with strong d^^'^-r^ character is smaller compared 
to dx^-y2. The wider (ia;2_j,2 -related bands are more oc- 
cupied than the narrower (i322_^2 band. The calculated 
eg orbital DOS for (LNO) i/ (LAO) i heterostructure is 
presented in Fig. [5l^a). Fig. (^fb), the integrated DOS of 
Fig-[2la), clearly shows that dx2_y2 orbital is more occu- 
pied than c?3j,2_r2. Note that the orbital polarization is 
present in the absence of strain. 

For a surface, the symmetry is also broken along the 
z-direction, but, in addition, the Ni coordination is re- 
duced. The calculated eg DOS of the single layer slab is 
shown in Fig. [UJc). First, the band width of the d^22_^2 
states is less than 2 eV which is not only significantly 
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FIG. 3: (Color online) Projected DOS for Ni Cg states in a 
LNO/LAO heterostructure ((a) and (b)), and single layer slab 
((c) and (d)). (a) and (c) are calculated under compressive 
strain (—3% in-plane lattice mismatch), and (b) and (d) under 
tensile strain (+3% in- plane lattice mismatch). Red (solid) 
and blue (dashed) lines correspond to d-j^2_j.2 and d^2_y2 
states, respectively. The Fermi level is set to be zero. 



smaller than the bulk value of ~ 4.5 eV, but also re- 
duced compared to the LNO/LAO heterostructure where 
the band width is larger than ~ 2 eV. On the other hand, 
the width of the (i3^2_j,2 -related bands remains compara- 
ble to the bulk value. For the heterostructure, the nar- 
rower d^z^-r^ band is located roughly in the middle of 
the dx2-y2 band, and the Fermi level is slightly below 
the middle point of the two bands, and the d^2_^y2 orbital 
is therefore more occupied than the d^^^^^^ (Fig. [SJa)). 
As a result, positive orbital polarization is induced as is 
clearly seen in Fig. [2jb). In the surface Ni site, on the 
other hand, the energy level of the d^z^_j.2 state comes 
down to the bottom of d^2_y2 band due to the change in 
nickel coordination and the Fermi level is roughly at the 
middle of d^^^-r^ state (Fig. [Sfc)). As a result, d^z^^r^ 
orbital is more occupied than dj.2_y2. This corresponds 
to a negative orbital polarization, see Fig. [^Jd). Again, 
this negative orbital polarization at the surface Ni state 
is a direct result of its local geometry and present without 
any additional strain effects. 

Let us now turn our attention to the effects of the 
strain. We took the relaxed LNO cubic bulk lattice pa- 
rameter (~ 3.81 A) as a reference, and applied compres- 
sive and tensile strain by setting the in-plane lattice pa- 
rameter to be smaller and larger, respectively. The ±3% 
of strain considered here is in a reasonable range to sim- 
ulate the experimental situation. For example, LAO, a 
widely used substrate for the compressive strain, has a 
lattice parameter of 3.78 A, and SrTiOa, a popular sub- 
strate material for compressive strain, has 3.91 A. 

Fig. |3l^a) and (b) shows Ni DOS in the heterostruc- 
ture under compressive (—3%) and tensile (-1-3%) strain, 
respectively. Similarly, Fig. (SJc) and (d) represents sur- 
face Ni DOS for —3% and +3% strain, respectively. The 
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FIG. 4: (Color online) Calculated orbital polarizations for a 
LNO/LAO heterostructure (blue), bulk LNO (green), and un- 
relaxed (red) and relaxed slab (orange) as a function of strain. 
Horizontal dotted line shows the zero polarization points. The 
relaxed result for heterostructure (magenta) is from Ref. [l^ 
with re-normalization. 



strain shifts the relative energy positions of the bands. 
However, the changes in orbital occupation are small 
and no sign reversal in the orbital polarization is found. 
The calculated polarization of the compressive and ten- 
sile strained heterostructure is 0.06 (Fig. |3l^a)) and 0.19 
(Fig-EIb)), respectively, and the polarization of surface 
is —0.25 (Fig. Eljc)) for compressive case and —0.12 for 
tensile (Fig. [Sjd)). 

The results of the changes in orbital polarization are 
summarized as a function of strain in Fig. |31 For bulk 
LNO, the two eg orbital states are degenerate under the 
cubic symmetry, and strain plays the major role in deter- 
mining the sign of orbital polarization. Under the —3% of 
compressive strain, c?3j,2_^2 orbital is more occupied lead- 
ing to negative orbital polarization, and the compressive 
strain results in the positively polarized orbital occupa- 
tion. Therefore the simple electrostatic picture works 
well in the bulk case. For heterostructures and surfaces, 
however, the symmetry breaking leads those systems to 
have asymmetric electronic structures and orbital occu- 
pations even without any strain as discussed in the above. 
The lattice strain is playing its role on top of the modi- 
fied electronic structures, and does not change the sign of 
polarization within the range of ±3% lattice mismatch. 
As is clearly seen in Fig. 21 in the heterostructure, Ni or- 
bitals are positively polarized and remain positive down 
to —3% strain. In the surface case, the negatively polar- 
ized orbital states is retained up to -1-3%. 

Rearrangements of atomic positions at the surface of- 
ten play an important role and cause a significant change 
in the electronic structure. It is therefore important 
to check the relaxation-induced change of the electronic 
structure and orbital occupation. From the force mini- 
mization calculations, it is found that the basic features 
of the electronic structure and orbital polarization at the 
surface do not change even after the relaxation. We per- 
formed the calculations within the three different in-plane 
lattice parameters corresponding to the compressive, zero 
and tensile strain, and the result of the calculated or- 
bital polarizations are presented in Fig.|3](dashed-dotted 
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FIG. 5: (Color online) Projected Ni-e^ DOS in the double 
layer slab structure, (a) represents to the top layer Ni site 
that has only five neighboring oxygens, and (b) shows the 
second layer Ni site that has six neighbors. Red (solid) and 
blue (dashed) lines correspond to d-j22_^2 and d^2_y2 states, 
respectively, and Fermi level is set to be zero. 



line with filled square boxes) from which one can see 
some deviation of the polarization trend but the ba- 
sic feature is retained. Ionic movements make polariza- 
tion slightly smaller and larger under the zero and -1-3% 
strain, respectively, with respect to the unrelaxed results 
(solid red line with open boxes). It is found that the 
relaxation-induced change becomes much larger in the 
compressively strained situation, which is attributed to 
the large ionic displacement along the out-of-plane direc- 
tion caused by the compressive in-plane strain, and the 
volume conservation often being observed in the related 
materials [H, [13 • ^^e zero strain case, the optimized 
distance between Ni and apical oxygen is ~ 2.04 A . 
It is decreased by 4% under the tensile strain and in- 
creased by 11% under the compressive strain. La ions 
also moved slightly in response to strains; toward and 
out-of the Ni02 plane by 0.04 and 0.09 A in case of ten- 
sile and compressive strain, respectively. In-plane oxy- 
gens (in the Ni02 plane) are found to move toward LaO 
plane by 0.06 A in case of compressive strain while it 
stays at almost same position under tensile strain. It 
should be noted that the calculated ionic movements is 
from the very thin slab calculation. In realistic situation 
in which the film thickness is much larger (> 10 nm) and 
it is on the substrate, the induced polarization change 
is much smaller and therefore is closer to the unrelaxed 
result (Fig. ID). 

Another important issue regarding the surface is the 
thickness of surface-like layers. To understand this point, 
we performed two-layer slab calculation. In this struc- 
ture, the first layer Ni ion {i.e., Ni(l) in Fig. [IJ has the 
reduced coordination, and the second layer Ni site {i.e., 



Ni(2) in Fig. [T]) has the octahedron oxygen environment 
as in the bulk. As shown in Fig. [5jb), the cubic sym- 
metry is significantly restored at the second-layer Ni site 
and the two eg states become bulk-like. The difference 
in the two eg orbital occupation is actually small, and 
the calculated orbital polarization is —0.02. The small 
deviation between (i3^2_^2 and dx2_y2 is the proximity 
to the surface state. The first layer Ni is similar to the 
single layer case (Fig. [Sja)), the polarization of which is 
—0.42. In addition, it is noted that the d^^2_^2 DOS at 
the Fermi level is reduced in Ni(2). Large DOS in Ni(l) as 
also observed in single-layer case may indicate the possi- 
ble unstability of the surface and the other reactions such 
as chemisorptions or stoichiometric changes. 

Our results shed new light on the interpretation of re- 
cent experiments. It is found from XLD measurement 
that compressively strained LNO thin film has negative 
orbital polarization whereas the tensile strain does not 
produce the sign reversal in the polarization [l7| . which 
is in a contradiction to the simple electrostatic Madelung 
picture assu ming the same energy levels both for dx2_y2 
and d 
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Our calculation provides a natural ex- 
planation as a surface effect dominant situation; i.e., the 
electronic structure change mainly caused by the reduced 
coordination is responsible for the negative polarization, 
and the effect persists even under the tensile strain. The 
case of heterostructure is more delicate especially be- 
cause the XLD shows no orbital polarization in the ten- 
sile strain case for which the simple electrostatic theory 
predicts the positive orbital polarization. According to 
our calculations, the tensile strain should enhance the 
positive orbital polarization in heterostructure. The ob- 
servation can be understood if there is some contribution 
from surface-like states even for the heterostructure case. 
As the dx^-y^ orbital is preferred in the heterostructure 
while d322_^2 is at the surface, there can be a compensa- 
tion between the two which may result in the negligible 
polarization along with the relatively small effect from 
the strains. 



IV. SUMMARY 

In conclusion, using density-functional theory, we have 
shown the change in the electronic structure of nicke- 
lates as a function of reduced dimensionality and strain. 
In surface Ni, the reduction of the coordination leads 
to a significant local crystal field. This leads to a nega- 
tive orbital polarization both in the absence and presence 
of strain effects. It is notably different from the posi- 
tive orbital polarization found in heterostructures where 
the asymmetric DOS caused by the strongly reduced hy- 
bridization along the out-of-plane direction results in a 
higher electron occupation of the x'^ — orbitals. Strain 
is found to play a relatively minor role on top of the mod- 
ified electronic structure of each system, and the sign of 
polarization does not change under the moderate strains. 
Our results are providing new insights for understanding 
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recent experiments on the related systems. 
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